First-principles study of orthorhombic CdTiOa perovskite 



(N 
O 

o 

(N 

> 
O 



o 



I 

C 

o 
o 



(N 

o 

O 

c3 



C 

O 
o 



13 



G. Fabriciusf] and A. Lopez Garcfajj] 
Depto. de Fisica, Univ. Nacional de La Plata, 
1900 La Plata, Argentina 
(Dated: February 1, 2008) 

In this work we perform an ab-initio study of CdTiOs perovskite in its orthorhombic phase using 
FLAPW method. Our calculations help to decide between the different cristallographic structures 
proposed for this perovskite from X-Ray measurements. We compute the electric field gradient 
tensor (EFG) at Cd site and obtain excellent agreement with available experimental information 
from a perturbed angular correlation (PAC) experiment. We study EFG under an isotropic change of 
volume and show that in this case the widely used " point charge model approximation" to determine 
EFG works quite well. 
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Perovskites materials of the form ABO 3 have been in- 
tensively studied last years. In addition to the interest 
that present the technological aplications of ferroelec- 
tric materials, perovskites are known to undergo several 
phase transitions as a function of temperature that make 
them appealing for theoreticians and experimentalists. 
The transitions between different structures usually in- 
volve very little atomic displacements and are difficult to 
characterize, requiring some times several experimental 
techniques. First principles calculations have been play- 
ing an increasing role in the understanding of the these 
systemsEJ. These studies not only try to give theoreti- 
cal explanations to the observed phenomena, but usually 
give valuable information that help in the interpretation 
of the experimental results. 

CdTiC>3 perovskite has been Jjhe. object of several 
studies using pfnRay diffractionHH'Q, dielectric prop- 
erties meihodsQtl'El, infrared spectroscopy^ and PAC 
tcchniquetm and there have been some difficulties and 
controversies in determining the symmetry group of 
it's orthorhombic structure at room temperature. In 
1957 H.F.Kay and J. L. Miles established through X- 
Ray measurements that CdTiOs perovskite should 
have oxthorhombic symmetry and non-centrosymmetric 
P6n2ili3 space group with a large displacement of Ti 
atoms along the b axiscl. However, the measurement 
of the electrical properties of CdTiOs single crystals 
couldn't resolve whether they were ferroelectric or not, 
then, the displacement of Ti atoms from their symmet- 
ric positions couldn't be confirmedQ. In Reference |8| 
Baumvol and coworkers study the EFG at Cd in CdTiOs 
through PAC measurements. Using the Point Charge 
Model approximation (PCM) to estimate theoretically 
the EFG they infer that the structure of the perovskite 
couldn't be the one corresponding to Pbn1\ space group 
but should be another one with centrosymmetric space 
group, also considered and discarded by Kay&Miles. 



Baumvol et al. made their arguments based on the asym- 
metry of EFG tensor, but in fact they obtained a dis- 
agreement of about 80% in the absolute value of EFG 
components. In 1987 S.Sasaki et alu performed new X- 
Ray studies on CdTiOs. They obtained that two or- 
thorhombic structures with very similar coordinates and 
different space groups (centrosymmetric Pbnm and non- 
ccntrosymmetric Pbn2\) fitted equally well the X-ray 
data. Analizing some structural and physical properties, 
they finally argued that it was more suitable to assign 
the space group Pbnm to CdTiOs. 

The main objective of the present work is to clarify 
this picture, studying from first-principles the electronic 
properties of the different structures proposed for the 
CdTi0 3 perovskite. 

Other important question that may be learned from 
this system concerns the correctness of the predictions 
made by FLAPW and PCM approximation to the EFG 
tensor at Cd. With PAC technique the EFG tensor at a 
probe introduced in the solid can be determinedtil, and 
Cd is one of the most used PAC probes, so, it is of valu- 
able interest to check the different theoretical approaches 
to the problem of determining EFG at a Cd site. In the 
PCM the EFG at the probe is usually approached as: 



V« = (l-7oo)V# tt - 
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where V^ 1 ' is the EFG produced by the valence charges 
located as point charges at the ions positions. In this 
approximation the only effect of the probe is to anti- 
shield the lattice contribution through the Sternheimer 
antishielding _£actoni3, that for Cd is usually taken as: 
7oo = — 29.27L3. PCM has been extensively use^,tar£¥ip-| 
diet EFG at Cd sites in perovskites and oxidescJtall3'cZl 
but calculations at the ab-initio level are rare in these 
systems because Cd usually enters the solid as impurity 
and this introduces several complications, as the effect of 
atomic relaxations, for examplola. In the case of CdTiOs 
perovskite, Cd is one of the constituents of the system, 
so, accurate ab-initio calculations can be performed and 
used to check the validity of PCM aproximation. 

The first-principles calculations were performed with 
the WIEN97 implementation of the Full-Potential 
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Linearized-Augmented-Plane=Wave (FLAPW) method 
developed by Blaha et al.E3. We use the General- 
ized Gradient Approximation (GGA) to describe the ex- 
change correlation potential in the parametrization of 
Perdew et altB. The size of the basis set in these 
calculations is controlled by the parameter RKmax 
that we fixed in the value of 8. The atomic sphere 
radii where chosen as 1.06, 0.90 and 0.85 A for Cd, 
Ti and O atoms respectively. Integrations in recipro- 
cal space were performed with the tetrahedron method 
using 18 and 36 fc-points in the irreducible hrst Bril- 
louin zone for Pbnm and Pbn2\ structures respectively. 
The unit cells have 20 atoms containing 5 molecules of 
CdTiOs. In the case of Pbnm structure there are 4 in- 
equivalent atoms in the cell: Cd(xcd,2/Cd, \), Ti(i,0,0), 
01(aroi,#oi,|) and 02(x O 2, 2/02,^02), while for Pbn2 1 
structure there are 5 inequivalent atoms: Cd(a;cd,2/Cd,:|), 
Ti(z T i,yTi,2Ti), 01(2;oi, 2/01, zoi), 02(a;o2, 2/02,^02) and 
03 (2:03, 2/03,-203)- Once self-consistency of the potential 
is achieved the Vij componentes of the EFG tensor are 
computed from the lattice harmonic expansion of the 
self-consistent potential^. The EFG tensor is then di- 
agonalized and it is expressed in terms of its greatest 
absolute value component: V33 and the asymmetry pa- 
rameter: 77 — (V 2 2 — Vn)/Va3, where Vu are the com- 
ponents of the diagonalized tensor that are ordered such 
that: \Vsa\ > \V 22 \ > \V U \. 

In Table | we show the results of our FLAPW calcu- 
lations performed for the different structures reported in 
the literature compared with the PAC determination of 
I V33 1 and 77. We see that the huge value obtained for the 
forces in the case of Pbn2\ structure of Ref. || and the 
poor agreement of eq and 77 with experiment lead us to 
discard this structure as the stable one for CdTiOs per- 
ovskite. In the case of structures of Ref. ||, we see that 
for both of them we obtained eq and 77 in excellent agre| 
ment with the experimental result of Baumvol et al. 
The forces obtained for these structures are much smaller 
than those obtained for Pbn1\ structure of Ref. || but the 
ones on oxygen atoms for Pbn1\ structure of Ref. || are 
still too big to be attributed to the precision of our cal- 
culations. So, in order to check the stability of structures 
of Ref. ^ in the frame of FLAPW calculations, we have 
relaxed both structures (keeping fixed the lattice param- 
eters at their experimental values) until forces on every 
atom are below a tolerance value taken as O.OleV/ A. 
These relaxation processes involve the variation of the 7 
and 14 independent coordinates that define Pbnm and 
Pbn1\ structures respectively, and we have performed 
them following a Newton damped schemec3. At the end 
of the relaxations we obtained that both structures re- 
laxed to the same Pbnm structure within a tolerance in 
distances of 0.01 A. The finally obtained Pbnm structure 
was very similar to the one proposed by Sasaki et al., as 
expected, due to the small forces previously obtained for 
this structure shown in Table |. So, our calculations in- 
dicate that Pbn1\ structure is not stable and confirm 
the assignation of Pbnm as the space group of CdTiOs 



perovskite that was suggested in Ref. |[ The fact that 
X-Ray data from Ref. [| are also adjusted quite well by a 
non-centrosymmetric structure with oxygens 02 and 03 
that slightly appart from their corresponding symmetric 
positions in Pbnm structure may indicate the presence 
of a soft mode in Pbnm structure involving those dis- 
placements. A future study of the phonons in the Pbnm 
structure would be valuable in order to check this hy- 
pothesis. 

As was mentioned, Baumvol et al. analizing the asym- 
metry of the EFG tensor with the point charge model 
supposed CdTiOs should have a centrosymmetric space 
group, but in fact they obtained V33 larger than exper- 
iment, at least, in 80 %. In effect, the discrepancy in 
V33 would be only about 80 % assuming the same sign 
for PCM prediction and experimental value, and they 
atributted this discrepancy to a failure of PCM to ac- 
count for a suppossed existing covalent contribution to 
EFG at Cd in this system. But the sign of V33 predicted 
by Baumvol is opposite to the one that emerges from 
FLAPW calculations. The point is that they used coor- 
dinates of a centrosymmetric structure studied in Ref. |^. 
If the coordinates obtained later by Sasaki et ala are 
used in the PCM calculations, V 33 =-4.41 and 77=0.45 
are obtained in fairly good agreement with experiment 
and FLAPW calculations. In order to inspect further 
this agreement between PCM and FLAPW predictions 
to EFG tensor, we perform FLAPW calculations for an 
isotropic change of volume of the Pbnm structure. For 
each volume considered we relaxed the internal coordi- 
nates and evaluate EFG tensor at the relaxed coordi- 
nates. The change is isotropic in the sense that we hold 
the relations a/b and a/c fixed in the experimental ra- 
tios of Ref. H The results are shown in Table [n]. Both 
EFG estimations predict an increase of V33 and 77 with 
compression and also an increase of the angle between 
V33-direction and x-axis. PCM not only predicts the 
same tendence as FLAPW, but also the absolute values 
are in fairly good agreement. However, the agreement of 
V33 and especially 77 seems to be better for larger Cd- 
O distances where covalency between Cd and oxygen is 
smaller. In Ref. |l5|, using PCM to compare with several 
PAC measurements in oxides, the authors obtained em- 
pirically that when Cd-0 bond is larger than 2.1 A a good 
agreement between predictions and experiment is found. 
The results of Table [il] suggest that the previous affirma- 
tion could be generalized including pcrovskites, since the 
agreement of PCM with FLAPW begins to worsen when 
d nn (Cd — O) approaches this value. 

In their PAC study of CdTiOs as a function of temper- 
ature, Baumvol et al. obtained that V33 and 77 increase 
when T decreases. This confirms the tendency predicted 
in our calculations assuming a positive thermal expansion 
coefficient for CdTiOs. The augment of 77 with pressure 
may be explained because the distorsion of the perovskite 
structure also increas es w ith pressure. This can be seen 
by looking at Tabic III where we show the behaviour 
of three parameters defined in Ref. to caracterize the 
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TABLE I: Absolute value of forces on inequivalent atoms (F) in eV/ A, component of largest absolute value of the EFG tensor 
(V33) in 10 21 V/m 2 and asymmetry parameter (77) calculated in the present work compared with PAC experiment. 







F(Cd) 


F(Ti) 


F(Ol) 


F(02) 


F(03) 


V33 


1 


Pbn2 1 : 


Ref. 


| 1.84 


3.43 


3.85 


1.04 


3.63 


+10.00 


0.80 


Pbnm: 


Ref. 


I 0.03 


0.00 


0.05 


0.05 




-5.10 


0.42 


Pbn2v. 


Ref. 


§ 0.03 


0.10 


0.41 


0.45 


0.54 


-5.14 


0.42 



PAC exp.B 5.41±0.80 0.407±0.008 



TABLE II: Electric field gradients as a function of an isotropic change of volume for the Pbnm structure of CdTiOs perovskite. 
PCM calculations were performed on the FLAPW-relaxed structures. V/Vo is the relation between volumes taken in our 
calculations and those of Ref. g d nn (Cd~0) (in A) is the nearest-neighbour Cd-O distance. V33 and 77, as in Table Q. 



FLAPW PCM 



V/Vo 


d nn (Cd-0) 


V33 


V 


V33-direction 


V33 


n 


V33-direction 


0.913 


2.18 


-5.63 


0.49 


(0.759, 0.651, 0.) 


-4.67 


0.59 


(0.773, 0.634, 0.) 


0.956 


2.21 


-5.33 


0.44 


(0.772, 0.636, 0.) 


-4.53 


0.52 


(0.778, 0.629, 0.) 


1.000 


2.25 


-5.04 


0.40 


(0.787, 0.617, 0.) 


-4.35 


0.47 


(0.784, 0.620, 0.) 


1.046 


2.28 


-4.82 


0.40 


(0.807, 0.591, 0.) 


-4.21 


0.45 


(0.795, 0.607, 0.) 


1.093 


2.31 


-4.62 


0.38 


(0.829, 0.560, 0.) 


-4.07 


0.41 


(0.809, 0.587, 0.) 



TABLE III: Parameters showing the distorsion in CdTiOs 
perovskite as a function of an isotropic change of volume. 
The observed tolerance factor (t ba) and the tilting angles of 
the octahedra in ab-plane (4> a b) and bc-plane (<j>bc) as defined 
in Ref. 0. t obs = (Cd- 0)/(Vt(Tj - O)), where (Cd-O) and 
(Ti — O) are the mean interatomic distances with twelve and 
six coordination for Cd and Ti sites respectively. 



V/Vo 



tobs 



4>ab 



4>bc 



0.913 
0.956 
1.000 
1.046 
1.093 



0.9759 
0.9765 
0.9772 
0.9777 
0.9783 



141.52 
142.18 
142.87 
143.53 
144.24 



105.20 
105.11 
104.97 
104.81 
104.66 



distorsion of this and other perovskites. For a cubic per- 
ovskite t oos — 1, <f>ab — 180° and 4>bc = 90°, and when 
larger is the separation from these values, larger is the 
distorsion of the structure and the tilting of the octahe- 
dra. We have obtained that the mentioned parameters 
vary almost linearly for an isotropic change of volume of 
CdTiOs perovskite. 

Even when in this study we keep fixed the a/b and 
a/c relations, and so we don't obtain an unconstrained 
energy vs volume dependence, it is interesting to look 
at this E(V) curve compared with the one of the cubic 
phase (see Figure 1). The high energy difference obtained 
between both phases (two orders of magnitude larger 



than the one obtained for the same phases in BaTiOsEj) 
may explain wiw the transition to cubic CdTiOs has not 
been reportecfa. In fact, perovskites that crystallize in 
an orthorhombic Pbnm structure at room temperature, 
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FIG. 1: Energy vs. Volume for cubic and orthorhombic phases 
of CdTiOs perovskite. Energy is referred to the minimum 
obtained for the orthorhombic phase, and volume is in terr 
of the experimental volume of the orthorhombic phase, Vd 
Arrows indicate the positions of the minima of the plotted 
curves. Full circles: calculated points, lines: cubic- fitt. 



as CaTiOs, SrZrC>3 and SrHf03 undergo three transi- 
tions when the temperature is increased: orthorhom- 
bic (Pbnm) — > orthorhombic (Cmcm) — > tetragonal 
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(1 4/ mem) — > cubi cBii. For the case of SrHf03 per- 
ovskite, for example, the transitions occur for temper- 
atures around 870 K, 1000 K and 1360 K. As far as 
we know, no structural study of CdTiOa at high tem- 
perature has been performed. Moreover in PAC stud- 
ies of orthorhombic CdTiC>3 up to 1270 Kino signature 
of any phase transition has been reporteda. In a previ- 
ous study on SrHfOsEa a difference of 0.27 eV has been 
obtained between orthorhombic (Pbnm) and cubic struc- 
tures. Thus looking at the 0.8 eV obtained in the present 
work for CdTiOa we would expect the transition to the 
cubic phase, if present, to occur at very high temperature. 
Future experimental work on CdTiC^ at high tempera- 
tures would be valuable to check this prediction. 

In summary our calculations help to decide between 
the proposed structures for orthorhombic CdTi03 per- 
osvskite. We have obtained that within the precision of 
our FLAPW calculations the proposed Pbnm structure 



from Ref. || is stable while Pbn2\ is not. We have ob- 
tained very good agreement between the calculated EFG 
and experiment and have shown that PCM approxima- 
tion works quite well in this case. However, care should 
be taken when applying PCM to other systems with lower 
Cd-0 bond, since the agreement seems to worsen in this 
direction. Our energy calculations suggest that the tran- 
sition to the cubic phase, if present, would take place at 
quite high temperature. 
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